function [ VX ] = Variance( X , w )

w = w / sum(w) ;

EX = X' * w ;

VX = ( X.^2 )' * w - EX^2 ;

end

